Finite Temperature Behavior of the v = 1 Quantum Hall Effect in Bilayer Electron 
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An effective field theoretic description of u = f bilayer electron systems stabilized by Coulomb 
repulsion in a single wide quantum well is examined using renormalization group techniques. The 
system is found to undergo a crossover from a low temperature strongly correlated quantum Hall 
state to a high temperature compressible state. This picture is used to account for the recent 
experimental observation of an anomalous transition in bilayer electron systems (T. S. Lay, et al. 
Phys. Rev. B 50, 17725 (1994)). An estimate for the crossover temperature is provided, and it is 
shown that its dependence on electron density is in reasonable agreement with the experiment. 
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I. INTRODUCTION AND SUMMARY 

The novel structure of bilayer electron systems 
(BLESs) has recently garnered considerable attention 
|0-^| . These systems have proven to be appropriate can- 
didates for probing quantum phases of electron systems, 
such as the quantum Hall effect (QHE). A BLES can 
be made in a double quantum well, where a high and 
hard-wall barrier leads to formation of two separate lay- 
ers of electrons. It may also be realized in a single wide 
quantum well (SWQW), in which the barrier separating 
the electron layers originates from Coulomb repulsion of 
electrons in the well 0]. 

The low energy effective theory for a BLES can be de- 
scribed as an XY spin model with an in-plane magnetic 
field, in which the spin-ordered phase corresponds to a 
correlated quantum Hall state ||,|7],||]. In this spin anal- 
ogy, local spin orientations encode the coherence between 
the two layers, the spin stiffness is a measure of loss of 
exchange and correlation energy corresponding to spatial 
variation of the relative coherence, and the magnetic field 
is given by the tunneling amplitude. In the absence of the 
tunneling term, the system exhibits a finite temperature 
Kosterlitz-Thouless (KT) phase transition |Hj, which de- 
stroys the phase coherence and thus the QHE p-ffl] . For 
any non-zero (but small) value of the tunneling param- 
eter, however, the transition is known to be smoothed 
out into a crossover, from a low temperature saturated 
ferromagnetic regime to a high temperature disordered 
regime g|0g|. 

Recently, Lay et al. have reported an experimental 
observation of a finite temperature quantum Hall "phase 
transition" in a BLES, realized in a SWQW @. They 
find that collapse of the QHE gap takes place at a tem- 
perature T* , whose value decreases when the electron 



areal density N s is increased. 

Here we study the effect of thermal fluctuations on 
the quantum Hall behavior of a BLES in a SWQW. We 
sketch the derivation of the effective field theory for this 
system and argue that an additional Ising-like symmetry- 
breaking term should be added to the Hamiltonian. We 
provide numerical estimates for the coupling constants 
of the theory using Hartree-Fock (HF) variational wave- 
function of the many-particle system combined with the 
self-consistent local spin density approximation (LSDA) 
for single particle states. We then study finite tempera- 
ture behavior of the system using renormalization group 
(RG) techniques. We argue that the system undergoes 
a crossover from a low temperature ordered regime to a 
high temperature disordered regime and that the experi- 
mentally observed A s -dependence of T* is related to the 
monotonically decreasing iso-spin stiffness with increas- 
ing electron density in the SWQW. 



II. EFFECTIVE HAMILTONIAN 

The low energy state of a BLES at v = 1 can be well 
described by the following variational wave function 
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where the up or down chiral iso-spin states denote local- 
ized electrons in the left or right side of the well, respec- 
tively, and A" is a quantum number such as the Landau 
gauge orbital guiding center. The phase angle field <j>(X) 
denotes the relative local coherence between the two elec- 
tron layers, and is well known to entail lowest energy ex- 
citations of the BLES [R,p|. In Eq.(l) we have neglected 



the real-spin degree of freedom of the BLES. This ap- 
proximation is justified by the self-consistent LSDA cal- 
culations yielding a unique, fully spin-polarized solution 
at v = 1 for the whole range of electron densities used 
in the experiment jlj. The calculated Zeeman splitting, 
A 2 , is strongly enhanced by the exchange-correlation en- 
ergy and always larger than the gap between two lowest 
energy levels in the SWQW, as shown in Table |. 

The general form of the effective Hamiltonian for a 
BLES in a SWQW has been obtained in Ref. ||], where 
it has been assumed that the BLES can be described 
effectively by the Hilbert space of the two lowest sub- 
bands. The effective energy functional can be written 
as 7i e g = Vhf + T, where T is the tunneling term 
and the microscopic HF Coulomb energy of the BLES 
Vhf = (*|^|*) is obtained as 
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with m(X) = (cos(j)(X), sm<j)(X)). The Coulomb energy 
has been evaluated using the Landau gauge, where 
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TABLE I. The SWQW parameters for various densities 
iV s (10 1:L cm -2 ) at v — 1, and given well width of 750 A corre- 
sponding to the sample used by Lay et al. |y, obtained from 
Eq.flq) and self-consistent LSDA one-particle states. 



SWQW. For the BLESs with large enough overlap inte- 
gral in the middle of the well ]/lTH ; anc [ thus K 2^ are 
significant. Values of the couplings p s , t and k 2 , of the 
layer separation d, and of the Zeeman splitting calcu- 
lated for the sample used by Lay et al. Q are presented 
in Table | 

The first term in Eq. (|5|) represents the usual, rotation- 
ally invariant superfluid exchange coupling, which yields 
a spontaneous phase coherent state (^,g| . (Note that the 
field <p is compact.) The above Hamiltonian is thus an 
XY-model with two symmetry-breaking terms: the sec- 
ond and third terms in Eq.(0) represent a uniform in- 
plane magnetic field and an Ising-like anisotropy, respec- 
tively. In the absence of the tunneling and the Ising-like 
terms, the system exhibits a finite temperature KT phase 
transition p0| , which destroys the phase coherence and 
thus the QHE §,f7|. 
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are the appropriate form factors, e is the dielectric con- 
stant of the host semiconductor, and £q = ^Jhc/eB is the 
magnetic length. 

In the continuum limit, the effective energy functional 
has a gradient expansion that reads 
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in which we have (only) neglected higher derivative 
terms. The origin of the iso-spin stiffness p s , and k 2 that 
corresponds to "pair hopping" , is the HF loss of exchange 
energy 
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and t is the tunneling amplitude, defined as one-half 
of the spacing between two lowest energy levels in the 



III. RENORMALIZATION GROUP APPROACH 

To study the XY model with the symmetry breaking 
terms, we follow closely the approach by Jose et al. [ pT| |. 
Finite temperature behavior of the system is described 
by the partition function 
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The essence of the symmetry-breaking terms can be cap- 
tured by performing Villain expansions, which introduce 
discrete Coulomb-like charge species, called type-1 and 
type-2 vortices [ [Lljjl2| ] . On the other hand, the compact- 
ness of the field 4> is accounted for by introducing yet an- 
other type of vortex, called type-0 vortex, which describes 
singular behavior of an otherwise noncompact field <j>'; a 
combination that serves as a substitute for <j> [ jL2]|l3| ] . The 
latter vortex introduces a new coupling constant, namely 
the fugacity yo — e ~ Ec / kBT , which is controlled by the 
core energy of a type-0 vortex E c = 2i:p s e c , where e c is 
a numerical constant close to 1 p4J. The corresponding 
three-species coupled Coulomb gas can then be studied 
using standard renormalization group (RG) techniques 




FIG. 1. RG flows in the different domains of the parameter 
space. 

pjfl . The RG flow equations for the dimensionless cou- 
pling constant K = p s /kBT, and the "fugacities" yo, 
Hi = t/8TT 2 p s , and 2/2 = k 2 /4:TT 2 p s are obtained as 
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to the leading order in the fugacities. A scaling argument 
shows that the coupling j/2 is always subleading compared 
to y\. We thus set it to zero for the moment to simplify 
the RG picture of the problem, and will comment on the 
effects due to a nonzero 2/2 later. The RG flow diagram 
for Eq.(||) is shown in Fig. [j]. 

The system described by Eq.(|5|) (for n 2 — 0) is known 
to have a duality in the parameter domain where the Vil- 
lain expansion is applicable. The duality maps the high 
temperature regions to the low temperature ones, with 
the roles of the fugacities being exchanged and temper- 
ature being inverted |O,0|. This necessitates the ex- 
istence of self-dual points, i.e. those which remain in- 
variant under the duality transformation, at some inter- 
mediate temperature. It is easy to see from Eq.(H) that 
K^ 1 = 2ir and j/o = yi define the line of self-dual points, 
which is denoted as path 1 in Fig. pi. We note that in 
Fig. [0, we have only sketched the RG flow structure for 
K^ 1 < 2ir, and that the corresponding behaviors in the 
K^ 1 > 2ir region can be understood using this duality. 

The above RG equations show that there are three 
different regions in the parameter space: (I) < K~ x < 
7r/2, in which yo is irrelevant and y\ is relevant, (II) 7r/2 < 
K^ 1 < 8ir, in which both yo and yi are relevant, and (III) 



K^ 1 > 8tt, in which yo is relevant and 1/1 is irrelevant. 
In region I, the symmetry breaking term is dominant, 
as illustrated by path 2 in Fig. [I], and the system is in 
a locked-in regime. While type-0 vortices are bound in 
pairs, type-1 vortices are unbound. On the other hand, 
region III (that is dual to region I) is a disordered regime, 
in which type-0 vortices are unbound and type-1 vortices 
are bound in pairs. 

Region II consists of four different areas. Path 3 shows 
a typical behavior for region Ha defined by 7r/2 < K^ 1 < 
2ir, and j/i > yo- While K~ x is monotonically decreas- 
ing under the RG flow, both yo & n d yi initially tend 
to increase with yi being always exponentially dominant 
compared to yo , until the flow eventually passes onto re- 
gion I. Although there are always some unbound type-0 
vortices, the overwhelmingly larger number of free type-1 
vortices ensures that the system is still in the locked-in 
regime. 

Region lib, which is defined by 7r/2 < K^ 1 < 2ir, and 
yi < yo, is a crossover domain. During the early stages of 
the RG flow, both y and y\ increase to approach the yo = 
yi surface, while K~ l is also increasing. The paths seem 
to be approaching the self-dual line asymptotically as if 
they are being attracted to an "intermediate temperature 
fixed point." However, at some point they will eventually 
exit the region, and find their ways to either the low 
temperature (such as path 4) or the high temperature 
(such as path 5) regions. The fact that the renormalized 
fugacities approach the yo = y\ surface means that there 
are about as many unbound type-0 and type-1 vortices, 
which explains why this is a crossover region. Finally, 
the behavior of regions lie (2ir < K^ 1 < 8tt, and y\ < 
yo) and lid (2ir < K^ 1 < 8tt, and yi > yo), can be 
understood using the duality. 

For an XY model in a magnetic field, a crossover from 
a low temperature locked-in regime to a high tempera- 
ture disordered regime has indeed been suggested in the 
literature by many authors gj§,|ll|,|l3l ■ However, it is not 
clear how this picture can be reconciled with the exper- 
imental observation by Lay et al., which is suggestive of 
a finite temperature phase transition in the same system 
|jj. To resolve this issue, we note that the presence of 
the self-dual line, and the fact that some RG flows in 
the crossover region are seemingly attracted to it, makes 
the crossover so rapid that it might "look like" a phase 
transition that is smeared out due to finite size effects. 

To make an estimate for the temperature at which this 
rapid crossover takes place, it is necessary to characterize 
the crossover more carefully |lq ] . Since in region II both 
fugacities are relevant, we need to worry about the valid- 
ity of the perturbative RG approach. Each RG equation 
describing exponential growth of a fugacity presumably 
breaks down at a length scale £ at which the fugacity 
is renormalized to unity. This correlation length simply 
characterizes the average separation between unbound 
vortices of the corresponding type, and is a measure of 



IV. DISCUSSION 
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FIG. 2. Comparison of the theoretical (solid line) and ex- 
perimental (triangles) crossover temperatures as a function of 
the BLES density. 



the number of free vortices [12]. Therefore, it is reason- 
able to assume that the crossover happens where the cor- 
relation length for type-0 vortices equals that of type-1 
vortices. 

Putting this together with the above picture for a rapid 
crossover, we assert that the crossover corresponds to 
an initial point, given by the bare parameters if _1 (0), 
yo(0), and ?/i(0), which is going to be renormalized to 
K- l (l*) ~ 2tt, and y (l*) ~ Vl (l*) ~ 1 via the RG flow, 
where I* = ln(£/a), and a is a microscopic length scale 
set by the core radius of a vortex. This criterion yields 
two equations for the two unknown parameters £ and 
the crossover temperature T* , which could in principle 
be solved numerically. 

Rather than elaborating on the numerical calculations 
that are not particularly illuminating, we attempt to 
make further simplifying approximations. We roughly es- 
timate the integrals as J dt{l/2 - irK{t)) ~ c l* (1/2 — 

ttK(0)) and /„** dt(l/2 - l/47rK(t)) ~ al*(l/2 - 
1/4ttK(0)), where the coefficients < cq < 1 and 
< ci < 1 describe the decay of the coupling constant 
K . Using this approximation, we obtain an estimate 
for the crossover temperature as 
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Values for T*, given by Eq.(H) with c = ci = 0.3 and 
e c = 1, are shown in Fig. |2j together with experimental 
data from Ref. JJJ . The magnitude of the crossover tem- 
perature is of the order of 1 K, and it decreases with in- 
creasing the BLES density, in reasonalbe agreement with 
the experiment. 



A nonzero value for k 2 does not change the RG picture, 
because the Ising-like term is always subdominant com- 
pared to the magnetic field, except for the low tempera- 
ture regime where the behavior of the system is governed 
by the ground state properties. While the Hamiltonian 
(ra) has only one minimum for t/n 2 > 4, a double-well 
structure with a local minimum at <j) = tt, a global min- 
imum at 4> = 0, and a maximum at <fi = cos~ 1 (— t/4n 2 ) 
develops for t/n 2 < 4. If the system is initially prepared 
to be at the metastable state, it will decay to the ground 
state with a rate ~ exp(—S[(j) c ]/h) where <j> c is the classi- 
cal solution, interpolating between the two minima [ JL7| . 
As one may see from Table pi the condition for the exis- 
tence of the metastable state (k 2 > i/4) is not satisfied 
for a typical BLES. However, it might be possible to re- 
duce the value of t by applying a tilted magnetic field 

Finally, we note that a full quantitative account of the 
experiment can not be achieved until one correctly takes 
into account all the details, possibly by performing an ex- 
act diagonalization calculation. However, we believe our 
effective field theory approach to the problem, provides 
an understanding of the nature of the observed "transi- 
tion" in the experiment, and is able to predict the correct 
order of magnitude and trend for the "transition" tem- 
perature. 
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